An adjustable Brownian heat engine 
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A microscopic heat engine is modeled as a Brownian particle in a sawtooh potential (with load) 
moving through a highly viscous medium driven by the thermal kick it gets from alternately placed 
hot and cold heat reservoirs. We found closed form expression for the current as a function of the 
parameters characterizing the model. Depending on the values these model parameters take, the 
engine is also found to function as a refrigerator. Expressions for the efficiency as well as for the 
refrigerator performance are also reported. Study of how these quantities depend on the model 
parameters enabled us in identifying the points in the parameter space where the engine performs 
either with maximum power or with optimized efficiency. The corresponding efficiencies of the 
engine are then compared with those of the endoreversible and Carnot engines. 
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We are in an era where devices including engines are 
being miniaturized to micron or even nano sizes. The 
physics of such microscopic engines is quite different from 
that of the macroscopic physics we daily experience. Due 
to their length scale, they are subjected to overdamped 
Brownian motion. Any task such an engine performs, 
such as translocating chemicals with finite velocity, costs 
energy. In order to find out how efficient such an engine 
is, one has to generalize the definition of efficiency. This 
is what Derenyi et al. did in their crucial paper on gen- 
eralized efficiency (l). Another issue of how an engine 
operates is whether a fastest transporting velocity exists 
under a given condition. Such fast transport is at the 
expense of losing significant energy and implying less effi- 
ciency. A compromise between fast transport and energy 
cost may lead to an optimized efficiency. Hernandez et 
al. recently came up with a unified optimization method 
for doing exactly this ||. 

The aim of this paper is to model an adjustable mi- 
croscopic (or Brownian) heat engine that could perform 
tasks in different modes of operation as mentioned above. 

The idea of a Brownian heat engine working due to 
nonuniform temperature first came up with the works of 
Biittiker, van Kampen, and Landauer Q while they were 
involved in exposing the significance of the now influen- 
tial papers of Landauer on blowtorch effect Q . Following 
Biittiker 's work, Matsuo and Sasa || considered such an 
engine as an example to show that it acts as a Carnot 
engine at quasistatic limit. Derenyi and Astumian 
after analyzing the details of heat flow in a Brownian 
heat engine, found that its efficiency can, in principle, 
approach that of a Carnot engine. The same group later 
used a model Brownian heat engine as an example to 
apply their newly introduced definition of generalized ef- 
ficiency and find its implications [jj] . 

All these recent works addressed how a Brownian heat 
engine behaves in the quasistatic limit. Instead of taking 
the engine as an example to address an issue of purely 



theoretical interest like quasistatic limit, we would like 
to take an exactly solvable model where information of 
practical interest can also be extracted from the model 
parameters themselves. Not only will we be dealing with 
the quasistatic limit but with general characteristics of 
real engines performing tasks either at maximum rate 
or at optimized efficiency or otherwise. Moreover, some 
of the wrong conclusions that were given due to either 
overlooking or lack of exact solution will be corrected. 

The model consists of a Brownian particle moving in a 
sawtooth potential with an external load where the vis- 
cous medium is alternately in contact with hot and cold 
heat reservoirs along the space (or reaction) coordinate. 
The shape of a single sawtooth potential, U s (x), located 
around x = is described by 



U s (x) = 



U 



1 



U + 1 



if -L x < x < 0; 
if < x < L 2 . 
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FIG. 1. Plot of the sawtooth potential in the presence 
of constant external load. The temperature profile is shown 
above the potential profile. 

The potential corresponding to the external load is lin- 
ear, fx, where / is the load. The temperature profile, 
T(x), within the same interval —L\ < x < L2 is describe 
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by 



T(x) = 



T h , if -Lx < x < 0; 
T c , ifO<x<L 2 . 



(2) 



Both U s (x) and T{x) are taken to have the same period 
such that U s (x + L) = U s (x) and T{x + L) = T{x) where 
L = L\ + L 2 . Note that the left side of each sawtooth 
from its barrier top overlaps with the hot region of the 
medium while the right side overlaps with the cold region. 
The sawtooth potential with the load, U(x) = U s (x)+fx, 
and the temperature profile, T(x), are shown in Fig. 1. 

From practical point of view, the size of such micro- 
scopic heat engine is limited by how small in size tem- 
perature gradients can be. Since temperature gradients 
of micron size can be produced [Q, the size of such an 
engine could be as small as few microns. 

Due to the presence of the hot and cold regions within 
each sawtooth potential and the external load, the Brow- 
nian particle will generally be driven unidirectionally and 
attain a steady state current, J, whose magnitude and 
direction depends upon the quantities characterizing the 
model. The dynamic equation governing the motion of 
the Brownian particle in such highly viscous and inho- 
mogeneous medium is the corresponding Smoluchowski 
equation first derived by van Kampen and later by Jayan- 
navar and Mahato fl§l: 



-i- (u\x)P(x,t) + |-(T(i)P(i,<)) 
J(X) \ ox 

(3) 



where P(x, t) is the probability density of finding the par- 
ticle at position x at time t, U'(x) = dU(x)/dx, and "f(x) 
is the coefficient of friction at position x. Boltzmann's 
constant, fee, is taken to be unity. The constant current 
at steady state is given by 
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U'(x)P ss (x) + — (T(x)P ss (x)) 



(4) 



where P ss (x) is the steady state probability density at 
position x. Using periodic boundary condition, P ss (x + 
L) = P ss (x), and taking coefficient of friction the same 
throughout the medium, we get a closed form expression 
for J for our potential and temperature profiles 



J = 



-F 



G\G 2 

where F, G\ and G2 are 
1, 

Gi = -4- (1 
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On the other hand, H can be put as a sum of three terms: 
A + B + C, where 
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* = ^(i-e-)(*-i), 
c = £(£)V-i-*>. 



(7) 



Note that a = (U + fL{)/T h and b = (U - fL 2 )/T c ) 
in Eqns. (6) and (7). For the case when there is no 
load (/ = 0) and L 2 = L\, the current takes a simple 
expression given by 



J = 



1 



2 7 (T h + T c 
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Note that this current is a resultant of currents to the 
right and to the left; i.e. J = J + — J_. 

One can clearly see that the model acts as a heat en- 
gine when the current is to the right. The particle un- 
dergoes cyclic motion wherein during each cycle it is first 
in contact with the hot region of width L\ and then with 
the cold region of width L 2 . We take account of energy 
flows between the two regions neglecting energy transfer 
via kinetic energy due to the particle's recrossing of the 
boundary between the regions PJlpC[|. As the particle 
moves through the hot region it absorbs a total amount 
of heat, Qh, which will enable it not only to climb up the 
potential (Uq + fLi) but also acquire additional energy 
■yvLi (v being the particle's drift velocity and equal to 
J{L\ + L 2 )) to overcome the region's viscous drag force 
so that 



Qh =U + ( 7 « + f)L x 



(9) 



On the other hand, the cold region will absorb energy as 
the particle moves down the potential hill while at the 
same time lose some energy due to the drag force in the 
region. Therefore, the net heat, Q c , absorbed by the cold 
region will be 



Qc = U - {iv + f)L 2 . 



(10) 



The net work, W, done by the engine in one cycle will 
then be the difference between Qh and Q c so that 



W=(jv + f)(L l + L 2 ). 



(11) 



Note that the term (7U + /) appearing in the above three 
expressions (Eqs. (9-11)) is the sum of the drag force due 
to the particle's motion through the medium and the load 
being lifted. One can take this term as generalized load so 
that the generalized efficiency, rj, as suggested by Derenyi 
et al. |l| will be given by 
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W_ = + f){L 1 +L 2 ) 

Q h Uo + irfv + f)^ ' 



(12) 



If, on the other hand, the load is large enough along with 
appropriately chosen other quantities then the current 
will be in the reverse direction in which case the load 
does work of amount Wl = f(Li + L 2 ) in one cycle forc- 
ing heat to be extracted from the cold region of amount 
Qc = Uq — (jv + f)L2- Under this condition the engine 
acts as a refrigerator. This will then lead us to a general- 
ized definition of coefficient of performance (COP) of the 
refrigerator, P re /, which is given by 



Pref — 



Qc _ Up - {jv + f)L 2 
W L /(Lr + L 2 ) 



(13) 



The condition at which the current changes its direction 
is the boundary demarcating the domain of operation of 
the engine as a refrigerator from that as a heat engine. 
In general, this condition is satisfied when 



U (T h -T c ) 
L\T C + L^Th 



(14) 



It is worth noting that the magnitude of the load at this 
point of zero current is exactly equal to what is usually 
called the stall force for molecular engines [[ll| . When we 
evaluate the expressions for both rj and Pref as we ap- 
proach this boundary, we analytically find that they are 
exactly equal to the values for efficiency of the Carnot en- 
gine and for COP of the Carnot refrigerator, respectively: 
T h -T c . ,. „ T, 



lim 77 



and lim P, 



T h ~ 1 -u ' ref 1,, /; 

this boundary at which current is zero corresponds to the 
quasistatic limit be it from the heat engine side or from 
the refrigerator side. We would like to emphasize that 
the efficiency of the heat engine at the stall force takes 
the maximum value contrary to the claim (of zero value) 
made by Derenyi, et al. We would also like to point 
out that the efficiency goes to zero as the barrier height, 
Uq, goes to infinity contrary to the claim (of maximum 
value) made by Derenyi and Astumian ||. This can be 
clearly seen from Eq. (12) above. We also claim that the 
expression for power of a reversible heat engine given in 
Eq. (10) of Derenyi et al. 0] does not make sense since 
one cannot have a thermal engine delivering finite power 
while, at the same time, operating reversibly. 



Hence, 
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FIG. 2. (a) Plot of j versus A for r = 1, 1=2 and u = 4. 
(b) Plot of j versus u for r = 1, 1=2 and A = 0.4. 



We next explore how current, efficiency and COP of 
the model depend on some of the quantities characteriz- 
ing it. In general, the quantities characterizing the model 
are Uq, Li, L 2 , f, T c and T h . We scale U , L 2 , T h and / 
such that u = U a /T c , £ = L 2 /L u r = (T h /T c ) - 1, and 
A = fL\/T c . Hence, we have four parameters u, I, r and 
A characterizing the model for a given T c and L\. We also 
scale current such that j = J/Jq, where Jo = T c /{~fL\). 
Figure 2(a) is a plot of the scaled current, j, versus scaled 
load, A. It shows that the engine works as a heat engine 
when the load value is less than the stall force (A = 0.8 
in this case) and works as a refrigerator when the load 
is larger than this value. One can also see the domains 
of operations of the model by plotting j as a function of 
u as has been shown in Fig. 2(b). Note that when the 
model works as a heat engine there is a finite u at which 
the current is maximum. This corresponds to the point 
in the parameter space at which the engine operates with 
maximum power. When we plot the efficiency, 77, versus 
A within the domain where the model works as a heat 
engine, we find that it increases linearly with increase in 
A attaining its maximum value (Carnot efficiency) at the 
stall force (See Fig. 3(a)). On the other hand, Fig. 3(b) 
shows plot of Pref versus A which starts from its maxi- 
mum value (COP of Carnot refrigerator) at the stall force 
and decreases fast as we increase the load. 
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FIG. 3. (a) Plot of n versus A for r = 1, 1=2 and u = 4. 
(b) Plot of P re f versus A for r = 1, 1=2 and u = 4. 

Let us now compare the efficiency of our model heat en- 
gine with that of the so-called endoreversible engine when 
both operate with maximum power. Curzon and Ahlborn 
|l2| took an endoreversible engine that exchanges heat 
linearly at finite rate with the two heat reservoirs and 
found its efficiency at maximum power, t/ca, to be equal 
to 1 — yj (Tc/Th). Figure 4 gives plots comparing t\ca 
with that of the efficiency of our model heat engine when 
it also operates with maximum power, T]mPi for different 
values of r. The plots show that the two results are rea- 
sonably close to each other for small r with rjcA being 
slightly smaller than t\mp- The two values coincide at 
a particular value of r and r\cA progressively gets larger 
and larger than that of t\mp for higher values of r. This 
shows that the assumption Curzon and Ahlborn took to 
get the expression for their efficiency holds reasonably 
well for our model heat engine as long as r is small. It 
is worth noting that our model did not take any extra 
assumption in the way it operates either with maximum 
power or otherwise. We, therefore, believe that the rea- 



3 



son for their deviation has got to do with the nonlinear 
nature of the way heat is exchanged. 
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FIG. 4. Plots of r/cA, tjmp, t/opt and tjcar versus r, 
where the model engine is put to function at / = and 1 = 2 
while u is fixed depending on whether it is working at either 
maximum power or optimized efficiency. 

We would further like to compare Carnot efficiency, 
Vcar, with what is called optimized efficiency, tjopt, by 
Hernandez et al. . This optimized efficiency is the effi- 
ciency at the point of operation of an engine where com- 
petition between energy cost and fast transport is com- 
promised. We briefly summarize the method Hernandez, 
et al. H used to get t)opt- When the model operates 
with finite time the amount of work, W, it delivers is 
such that it lies between the maximum, W max , and min- 
imum, Wmin , amount of work that can be extracted from 
the model: Wmm < W < W ma x- They defined two quan- 
tities: effective work, W e = W — W m i n and noneffective 
work W ne = W m ax — W and introduced a function, Q, to 
be optimized to be equal to the difference between these 
quantities. For our model heat engine (since W m i n = 0) 
this function is given by 
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FIG. 5. Plot of Q versus u for / = 0, I = 2 and r = 1. 

Figure 5 shows plot of Q versus u in which the function 
has indeed an optimum at a finite value of u. The effi- 
ciency of the model heat engine when it operates at this 
particular point in the parameter space is what we call 



as optimized efficiency, tjopt- For the sake of compari- 
son, we plot tjopt along with the corresponding Carnot 
efficiency, tjcar-i m the same figure (Fig. 4) as we did 
for the other efficiencies. The plots clearly show that 
rjopT always lies between the rjcAB. and T]mp- One can 
therefore say that the operation of the model at optimized 
efficiency is indeed a compromise between fast transport 
and energy cost. 

In conclusion, we believe that our work has, for the first 
time, quantitatively explored the basic properties such as 
current, efficiency and COP of a model Brownian heat 
engine. The recently introduced definition of generalized 
efficiency ||] has been used in getting the expressions for 
these quantities. We have found operating conditions of 
the model under which either maximum power is deliv- 
ered or optimized efficiency is attained and have been 
able to explicitly compare their efficiencies. These theo- 
retical results propose the extents and limits to be con- 
sidered in the design of actual Brownian heat engines 
which, of course, take finite time to accomplish any task. 
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